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1 INTRODUCTION 

The  laser  anemometer  is  an  instrument  for  measuring  the  dynamic  properties 
of  fluid  motion.  It  is  based  on  the  Doppler  effect  at  optical  frequencies  and 
uses  the  information  contained  in  the  temporal  variations  of  the  light  scattered 
by  particles  embedded  in  the  fluid.  It  offers  the  fluid  dynamicist  the 
capability  of  making  essentially  non-intrusive  studies  of  flow  situations  which 
were  previously  either  very  difficult  or  impossible  to  carry  out*. 

Early  laser  anemometry  measurements  at  RAE  were  made  in  naturally  seeded 
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flows  . However,  since  naturally  occurring  particles  can  vary  widely  m their 
characteristics,  including  size,  it  is  desirable  to  start  with  an  unseeded  flow 
and  to  add  to  it  particles  with  known  characteristics. 

In  the  closed  circuit  wind  tunnels  at  RAE  in  which  experiments  have  been 
carried  out  it  has  been  found  that  the  air  becomes  relatively  clean  after  a short 
running  time.  Consequently  it  only  remains  to  choose  a suitable  material  and  an 
appropriate  method  of  generating  seeding  particles  from  it.  There  are  several 
important  factors  to  be  considered  in  making  this  choice: 

(i)  If  the  particles  are  to  convey  meaningful  velocity  information  about  the 

flow,  they  must  faithfully  reproduce  the  motion  of  the  fluid  element  in 

which  they  are  embedded,  even  when  this  experiences  strong  forces;  for 

example  when  traversing  a shock  wave  or  following  a turbulent  wake.  There 

is  a large  amount  of  literature  on  this  subject,  but  there  are  still  some 

topics  which  are  only  partially  understood;  for  example,  particle  behaviour 
. . 3 

in  a rotating  and/or  turbulent  flow.  Boothroyd  gives  a general  review 
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of  the  subject  whilst  Karchmer  and  Kirsch  and  Mazumder  deal  with  oscil- 
lating flows.  Some  workers  use  the  behaviour  of  the  particle  in  strongly 
accelerating  flows  to  obtain  an  estimate  of  its  size;  see  for  example 
Refs  6-8.  The  conclusions  reached  in  these  studies  indicate  that  the 
particles  should  be  made  as  small  as  possible,  consistent  with  the  require- 
ment that  they  do  not  exhibit  any  Brownian  motion  in  response  to  the  motion 
of  the  gas  molecules. 

(ii)  The  light  backscattered  from  the  particle  must  be  sufficiently  strong  to 
make  experiments  possible  in  the  larger  wind  tunnels.  Clearly  for  the 
particle  to  scatter  any  light,  its  refractive  index  must  be  different  from 
that  of  the  host  fluid.  The  size  of  the  particle  also  affects  the 
scattering  efficiency. 
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(iii)  Toxicity,  abrasiveness,  cost  and  ease  of  generation  should  also  be  taken 

into  account. 

Oil  droplets  with  diameters  in  the  range  0. 1-1.0  micron  meet  these  cri- 
teria for  most  practical  applications  of  interest  to  RAE.  The  seeding  method 
adopted  during  the  laser  anemometry  development  programme  at  RAE  is  based  on  a 
commercial  oil-mist  generator,  designed  to  produce  a fine  oil  mist  suitable  for 
use  in  lubrication  systems.  Oil  consumption  is  found  to  be,  typically,  approxi- 
mately 2 ml  per  hour.  No  serious  problems  arising  from  wind-tunnel  contamination 
have  been  encountered. 

A variety  of  experiments  have  now  been  carried  out  using  this  form  of 
seeding;  of  particular  interest  to  studies  of  particle  behaviour  in  rapidly 
accelerating  flows  and  to  the  overall  accuracy  attainable  were  measurements  made 

, . . 9 

in  supersonic  flow  over  a ten-degree  cone  at  normal  incidence  . The  results  of 
this  experiment  could  be  compared  with  theoretical  calculations  of  the  velocity 
distribution,  for  a given  stagnation  temperature  and  Mach  number.  Experience  has 
also  been  gained  in  the  3ft  * 3ft  supersonic  wind  tunnel  at  RAE,  Bedford;  studies 
of  the  interaction  of  a normal  shock  wave  with  a turbulent  boundary  layer  gave 
quantitative  results'*^  which  accorded  with  theoretical  expectations  and  in 
regions  where  no  fully  detailed  theoretical  prediction  was  available,  new  informa- 
tion was  obtained. 

Although  these  experiments  have  verified  the  predicted  particle  behaviour 
an  independent  laboratory  study  of  some  of  the  particle  characteristics  was 
carried  out.  The  main  quantities  of  interest  were  the  refractive  index  and  the 
mean  diameter  of  the  droplets  produced  by  the  oil-mist  generator.  A specifica- 
tion of  the  oil  is  given  in  Ref  1 1 . 

A widely  used  method  for  the  determination  of  the  size  of  a small, 
spherical,  dielectric  particle  is  based  on  measurements  of  the  light  scattered 
by  the  particle  out  of  an  incident  beam.  Experiments  were  carried  out  to  deter- 
mine the  refractive  index  of  the  oil  and  the  angular  distribution  of  the  scattered 
light  field  for  a particular  wavelength  and  direction  of  polarisation  of  the 
incident  light;  the  effect  of  a settling  chamber  was  also  studied.  Finally  an 
attempt  was  made  to  verify  the  value  of  the  size  derived  from  the  results  of  the 
angular  scattering  experiment  by  measurement  of  the  turbidity  of  the  oil  mist  at 
several  different  wavelengths. 

This  Report  outlines  the  basic  theory  of  light  scattering  applied  to  small, 
spherical  particles  and  describes  the  experiments  carried  out  on  the  oil  droplets. 
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2 THEORY  OF  LIGHT  SCATTERING 
2. 1 Mie  theory 

Light  scattering  experiments  can  be  used  to  characterise  small  dielectric 

spheres  since  the  measurable  quantities  are  connected  by  exact  theoretical 

relationships.  The  original  formulation  of  these  relationships  is  due  to  Lorenz'' 

. 13 

although  they  are  normally  attributed  to  Mie  , who  discovered  them  independently 
some  years  later. 

Mie  theory  predicts  the  precise  form  of  the  electromagnetic  field  which 

results  from  the  scattering  of  a linearly  polarised  plane  wave  by  a spherical 
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particle  of  arbitrary  size.  The  theory  (see  for  example  Van  de  Hulst  or 
Kerker'^)  is  based  on  a solution,  subject  to  the  appropriate  boundary  conditions, 
of  a scalar  form  of  the  vector  wave  equation  arising  from  Maxwell's  electromag- 
netic theory.  The  origin  of  the  co-ordinate  system  is  taken  to  be  the  centre  of 
the  particle.  The  Mie  solution  enables  an  expression  to  be  written  which  des- 
cribes the  form  of  the  scattered  wave  in  terms  of  the  amplitude  functions  Sj(9), 
arising  from  the  components  of  polarisation  of  the  incident  wave  perpen- 
dicular and  parallel  respectively  to  the  scattering  plane: 


S,  (0) 


I 

n = 1 


2n  + 1 
n(n  + 1) 
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n sin  0 


P 1 (cos  0)  + b 
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„3?(pi<cos  e))} 
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t. 


The  scattering  plane  is  defined  as  the  plane  containing  the  incident  and  scattered 

wave  vectors,  and  the  scattering  angle  0 is  the  angle  shown  in  Fig  1.  The 

quantities  a and  b are  the  Mie  coefficients  and  the  P*  (cos  0)  are 
n n n 

associated  Legendre  polynomials.  The  Mie  coefficients  are  complex  quantities 
which  depend  on  the  particle  size,  refractive  index  and  the  wavelength  of  the 
incident  light. 

The  scattered  fields  E^  and  E^  , perpendicular  and  parallel  respectively 
to  the  plane  of  scattering,  are  given  at  a point  distant  r from  the  origin  in 
terms  of  the  incident  fields  E^  and  by: 
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E]  = S ] (0)  exp(ikz  - ikr)E|0)^ikr 


E2  = S2^6^  exP(ikz  " ikr)E2 


(0)/ikr  . 


(3) 

(4) 


Here  z is  the  direction  of  propagation  of  the  incident  plane  wave  and 
k = 27r/A  , where  A is  the  wavelength  of  the  light.  The  quantity  actually 
measured  in  any  light  scattering  experiment  is  the  intensity;  for  perpendicular 
polarisation 

1,(6)  = |E]|2  = ix(Q)lJk2r2 


where 


i,  (e)  = I s J (e) 


(5) 


and  for  parallel  polarisation 


I2(0)  = |E2|2  = i2(6)IQ/k2r2 


where 


i2(0)  = |S2(0) 


(6) 


For  naturally  polarised  light 


1(0)  = (i j (6)  + i2(0))IQ/2k2r2 


(7) 


1^  is  the  intensity  of  the  incident  wave. 

At  this  point  we  introduce  the  concepts  of  scattering  cross-section  and 

efficiency  factor.  Imagine  a plane  wave  incident  from  infinity  on  a spherical 

particle.  The  energy  scattered  in  all  directions  is  equal  to  the  intensity  of 

light  from  the  undisturbed  wave  which  would  have  fallen  on  some  area  C ; 

sea 

similarly  the  energy  absorbed  by  the  particle  is  equal  to  the  intensity  of  the 

undisturbed  incident  wave  which  would  have  fallen  on  some  area  . Finally 

let  the  total  energy  removed  from  the  incident  wave  be  equal  to  the  energy  of 

the  undisturbed  beam  over  some  area  C ; then  the  law  of  conservation  of 

ext 

energy  requires  that  the  following  relationship  is  true: 


c = c + C , 

ext  sea  abs 


(8) 


The  quantities  C , C and  C , have  dimensions  of  area  and  are  called 

cXl  SC3  cl  D s 

the  extinction,  scattering  and  absorption  cross-sections  respectively.  These 
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quantities  may  be  divided  by  the  cross-sectional  area  of  the  particle,  ira  , 

where  a is  the  particle  radius,  to  provide  the  dimensionless  quantities 

Q , Q and  Q , , known  as  the  extinction,  scattering  and  absorption 

ext  sea  abs 

efficiencies  respectively. 

The  magnitude  of  the  extinction  efficiency  can  be  determined  from  the  ampli- 
tude functions,  by  setting  the  scattering  angle  to  zero.  In  this  case 


S,(0)  = S.(0)  = S(0)  = 7 ) (2n+l)(a  + b ) . 

1 2 ^ n n 

n=  1 


It  can  be  shown  that 
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or 


Q = (4/a2)Re{S(0)}  , 

ext 


ext 


(2/a2) 


oo 

(2n  + 1 )Re{an 

n=  1 


+ b } . 
n 


(9) 


(10) 


a = 2na/X  is  the  particle  size  parameter.  Note  that  the  extinction  efficiency 
is  independent  of  the  direction  of  polarisation  of  the  incident  wave.  For  a non- 
absorbing dielectric  sphere,  with  an  entirely  real  refractive  index,  Q = 0 

and  the  scattering  efficiency  Q is  equal  to  the  extinction  efficiency  Q . 

sea  6xt 

It  will  be  shown  that  this  is  effectively  the  case  for  particles  of  interest 
to  this  study. 

The  forms  of  the  intensity  functions  and  the  efficiency  factors  are  extremely 

complex  and  require  the  facilities  of  a large,  fast  digital  computer  for  their 

evaluation.  However,  tables  exist  for  certain  combinations  of  size,  refractive 

index,  wavelength  and  incident  polarisation  direction.  This  Report  makes  use 

of  tables  of  Q (a)  , due  to  Penndorf*^,  calculated  for  spherical  particles 
sea 

with  a purely  real  refractive  index  equal  to  1.486.  It  will  be  shown  that  this 
value  corresponds  closely  to  the  actual  value  for  the  oil  used  in  the  experiments. 

A plot  of  this  curve  is  shown  in  Fig  2,  together  with  a similar  plot  for  a 
refractive  index  of  1.44;  this  shows  the  effect  of  a small  change  in  refractive 
index  on  the  form  of  the  curve.  The  intensity  functions  I ^ ( 0 ) , 12(6)  and 

I (0)  + I (0)  shown  in  Figs  3-6  were  derived  from  a computer  program  made 

. 2 1 7 . . . 

available  by  Dave  and  represent  typical  functions  for  a given  particle 

size  and  refractive  index  m for  a wavelength  X . If,  as  is  normally 
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the  case,  the  experiment  measures  the  composite  scattering  function  arising 
from  a collection  of  particles  with  differing  radii  then  the  finer  detail 
in  the  curves  for  a mono-dispersion  is  lost.  Under  these  circumstances  a 
practical  unambiguous  method  of  determining  the  particle  size  distribution  from 
the  experimental  data  is  not  known. 


2.2  Theory  of  the  angular  scattering  experiment 


One  method  of  overcoming  the  interpretation  problem  encountered  above  was 

18  2 
proposed  by  Sloan  who  showed  that,  if  the  curve  log  1(0)0  is  plotted  against 

the  logarithm  of  scattering  angle  0 , the  position  of  the  first  maximum,  0 * 

is  inversely  proportional  to  the  mean  of  some  linear  characteristic  dimension  of 

the  particle.  The  characteristic  dimension,  in  the  case  of  a spherical  scattering 

particle,  is  the  radius  a and  the  expression  which  relates  it  to  the  value  of 

0 is 

max 


a0  “9.2  micron  degrees 
max 


(11) 


The  constant  of  porportionality  (9.2  pm  degrees)  remains  approximately  valid  over 

the  range  of  size  from  0.2  to  100  microns.  Sloan  also  found  that,  to  a first 

approximation,  the  value  of  0 for  a particular  size  was  independent  of  the 

max  I g 

refractive  index  of  the  particle.  Mie's  theory  predicts  a departure  from 

equation  (11)  as  the  size  falls  below  10  microns  which  depends  on  the  refractive 

index.  A graph  showing  the  relationship  between  the  exact  and  approximate 
20 

values  is  shown  in  Fig  7. 

2.3  Theory  of  the  extinction  experiment 

Measurement  of  the  attenuation  of  a collimated  light  beam  at  different 

wavelengths  passing  through  a suspension  of  droplets  yields  information  about  the 
. 14 

particle  size  . The  variation  of  intensity  of  a beam  of  light  passing  through 
a cloud  of  uniformly  distributed  droplets  is  proportional  to  its  intensity  f 
so  that  the  extinction  occurring  in  the  distance  Az  in  the  direction  of  pro- 
pagation is: 

Af  = - xf Az  . (12) 


The  constant  of  proportionality  x is  the  extinction  coefficient.  Equation  (12) 
can  be  integrated  to  give 


f/f 


0 
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exp(-  xz) 


(13) 


where  is  the  intensity  of  the  light  at  the  point  from  which  z is  measured. 

This  equation  is  the  Bouger  or  Lambert-Beer  transmittance  law.  t can  be 

expressed  in  terms  of  the  extinction  efficiency  (Q  ) and  the  number  of  particles 

ext 

per  unit  volume  N : 

x = Nira^Q  . (14) 

ext 

If  the  suspension  of  droplets  comprises  a size  distribution  characterised  by  the 
function  p(a)  , so  that  p(a)6a  is  the  number  of  particles  per  unit  volume  with 
radii  in  the  range  a to  a + $a  , then  equation  (14)  becomes: 
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IT 
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2 

a p 


(a)Q  da 
ext 


(15) 


The  quantity  f/f^  can  be  measured  at  a number  of  different  wavelengths. 

Thus  if  the  form  of  the  function  p(a)  is  known  or  if  the  suspension  is  assumed 

to  be  mono-disperse,  the  curve  of  log^f^/f  plotted  against  the  wavenumber  v , 

defined  by  v = tt/ A , can  be  fitted  to  a similar,  calculated  plot  of  Q (a)  ; 

ext 

we  recall  that  a = 2ira/A  . Comparison  of  the  abscissae  yields  the  particle 

size,  which  in  the  case  of  a mono-dispersion  is  the  common  diameter.  In  the  case 
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of  a particle  size  distribution  it  is  the  Sauter  mean  diameter,  defined  as  the 
ratio  of  the  third  to  the  second  moments  of  the  size  distribution. 

3 EXPERIMENTS 

The  interpretation  of  the  results  obtained  from  light  scattering  experiments, 
described  in  this  section,  requires  a knowledge  of  the  complex  refractive  index 
of  the  oil. 

3. 1 Measurement  of  the  refractive  index  of  the  oil 

The  complex  refractive  index  m is  conventionally  written  in  the  form: 


m = n - in. 

r l 


(16) 


The  real  part  nf  is  related  to  the  velocity  of  propagation  of  an  electromagnetic 
wave  through  the  medium  and  is  the  part  used  for  the  calculation  of,  for  example, 
the  angles  of  refraction.  The  imaginary  part  n^  describes  the  attenuation  of 
the  wave  as  it  passes  through  the  medium  (ie  the  process  of  absorption) . The 
equation  for  the  electric  field  component  of  an  electromagnetic  wave  propagating 
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in  the  z direction  with  an  angular  frequency  w can  be  written  in  terms  of 
the  refractive  index  m : 

E = Eq  exp{-  iw(t  + mz)}  I 

E is  the  amplitude  of  the  electric  field  at  the  point  from  which  z is 
^ 2 

measured.  The  intensity  |e|  at  any  point  in  the  medium  is  given  by  the 


formula 


IQ  exp (-  Kz) 


where  1^  is  the  quantity  |Eq|  and  K is  the  absorption  coefficient: 


K = 4irn^/X^  , 


Aq  is  the  wavelength  of  the  radiation  in  vacuo. 

Measurement  of  the  quantity  I/I^  at  a given  wavelength  A and  over  a 
distance  z yields  the  imaginary  part  of  the  refractive  index  n^  . The 
schematic  layout  of  such  an  experiment  is  shown  in  Fig  8.  It  is  important  to 
restrict  the  power  density  in  the  test  cell,  otherwise  distortion  of  the  laser 
beam  may  occur  through  absorption  of  energy  and  consequent  heating  of  the  oil. 

The  distortion  will  arise  from  the  local  change  in  the  real  part  of  the  refractive 

index  n . Note  that  this  quantity,  for  liquids,  is  very  sensitive  to  tempera- 

r -4  . 22 

ture;  for  example  An^/AT  for  water  is  —10  per  degree  centrigrade  . 

The  value  of  the  real  part  of  the  refractive  index  was  determined  from  an 
experiment  with  a hollow  prism  in  a simple  spectrometer.  The  prism,  which  con- 
tains the  oil,  was  constructed  from  microscope  slides;  the  geometry  is  shown  in 

• A # A # 

Fig  9.  Measurement  of  the  angles  of  incidence  a and  refraction  8 would  yield, 

A 

via  Snell's  law,  the  quantity  n^_  if  direct  measurement  of  the  angle  B were 
practicable.  Using  the  optical  arrangement  shown  in  Fig  10  it  is  possible  to 
measure  the  angles  2a  and  p and,  with  an  accurate  knowledge  of  the  internal 

A 

angle  B , it  is  possible  to  evaluate  the  refractive  index  from  the  following 
relationships: 

A .A  A a A 

n^_  cos  B = (sin  5 + cos  B sin  a)/sin  B (20) 

A A 

sin  B = (sin  a)/n  (21) 


B - a + p - it 


L 


6 


(22) 


1 1 

A # A A 

The  internal  angle  B was  first  determined  by  measuring  the  angles  2a  and  p 
. 22 

for  a liquid  with  a known  refractive  index  (distilled  water)  and  then  using 
relationships  (20)  to  (22)  for  this  value  of  n^  . 

In  practice,  because  the  oil  and  water  are  immiscible,  both  experiments 
could  be  carried  out  together  and  a set  of  angles  p obtained  for  a constant 
angle  of  incidence  a and  a selection  of  wavelengths  (577,  546  and  436  nm) , for 
both  the  oil  and  the  water. 

Although  subsequent  light  scattering  experiments  employed  lasers  exclusively, 
a mercury  vapour  source  was  used  for  these  measurements,  mainly  to  avoid  the 
diffraction  effects  associated  with  highly  coherent  light. 

The  results  for  the  real  part  of  the  refractive  index  were  extended  to 

the  required  laser  wavelengths  632.8,  514.5,  488.0  and  457  nm  by  use  of  the 
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Cauchy  dispersion  relationship  : 

B 

n = A + — =■  + 

r A2 

3.2  Angular  intensity  measurements 

The  Mie  theory,  outlined  above,  gives  the  form  of  the  intensity  functions 
ijO),  ^(S)  which  would  arise  from  scattering  by  a single  particle.  The  deter- 
mination of  these  quantities  would  require  the  selection  and  positioning  of  the 

particle  while  measurements  were  carried  out.  One  possible  method  of  achieving 

24 

this  would  be  an  adaptation  of  the  Millikan  oil  drop  experiment,  although  this 
would  require  a complex  position-control  system  and  an  optical  enclosure,  which 
could  give  rise  to  parasitic  radiation. 

An  alternative  technique  is  described  in  Ref  25  which  employs  a high  speed 
rotating  mirror  to  change  the  scattering  angle  6 and  accumulates  data  from  the 
light  scattered  by  particles  which  pass  slowly  through  the  sample  volume.  This1 
system  works  over  a limited  range  of  6 and  also  requires  the  integration  of 
the  scattered  light  from  many  successive  particle  transits  to  provide  adequate 
signals  for  analysis. 

The  method  adopted  here  avoids  the  use  of  an  enclosure  and  therefore  the 
problems  associated  with  parasitic  radiation,  by  analysing,  with  a collimated 
detection  system,  the  light  scattered  by  a fine  stream  of  particles,  illuminated 
by  a focused  laser  beam.  A schematic  diagram  of  the  experiment  is  shown  in 
Fig  II.  Note  that  in  a typical  laser  anemometry  application  the  seeding  particles 


055 


12 


are  introduced  into  the  gas  in  a similar  manner,  and  also  that  during  the  course 
of  an  anemometry  experiment  the  detector  is  exposed  to  the  scattered  light  field 
arising  from  the  whole  particle  size  distribution. 

The  basic  optical  system  consists  of  three  parts,  an  illumination  system, 

a detector  arm  that  can  move  in  the  scattering  plane  over  a range  of  6 and  a 

static  reference  arm.  The  first  part  uses  the  radiation  from  an  argon-ion  laser 

(Spectra-Physics  Model  165),  operated  in  the  constant  power  output  mode  at 

approximately  200  mW.  The  output  beam  is  expanded  by  a lens  and  focused  to 

a waist  by  lens  L 2 . The  diameter  of  the  waist,  approximately  10  microns,  was 

calculated  from  measured  input  beam  characteristics  using  the  expressions  derived 
26 

by  Kogelnik  and  Li  . Measurements  of  the  diameter  of  the  laser  beam  as  it 
expanded  away  from  the  waist  confirmed  this  value.  Plots  of  the  beam  profile 
and  the  axial  power  distribution  are  shown  in  Fig  12.  The  output  of  the  laser  is 
vertically  polarised  by  virtue  of  the  orientation  of  the  Brewster  windows  within 
the  cavity.  Hence,  in  order  to  measure  I (0)  with  the  same  experimental 
arrangement  a rotation  of  90  degrees  of  the  direction  of  polarisation  direction 
is  required.  This  was  achieved  using  a half  wave  plate  Pj  . 

The  second  part  of  the  optical  system,  the  movable  detector  arm,  consisted 
of  a long  collimating  tube  which  could  describe  an  arc  in  the  horizontal 
scattering  plane,  centred  on  the  intersection  of  the  axes  of  the  particle  jet 
and  the  illuminating  laser  beam.  The  range  of  6 covered  was  from  10  to 
160  degrees.  It  was  not  possible  to  reduce  the  lower  limit  because  the  rapidly 
expanding  laser  beam  beyond  the  sample  volume  interfered  with  the  scattered  light 
signal.  The  value  of  the  scattering  angle  was  determined  from  a large,  graduated 
rotating  table  which  formed  part  of  the  mechanical  arrangement. 

The  collimating  tube  defined  a collection  angle  of  approximately  6.7  mrad 
(^.38  degrees)  by  means  of  two  1mm  diameter  apertures  separated  by  a light  tight 
tube  approximately  300  mm  long.  Attached  to  the  rear  of  this  tube  was  a photon 
counting  photomultiplier  tube,  type  ITT  FW  130,  fitted  with  a narrow  band  filter, 
matched  to  the  laser  wavelength  (488  nm) , and  an  imaging  lens  which  focused  the 
rear  aperture  of  the  collimating  tube  on  to  the  photo-cathode.  The  photo- 
multiplier tube  is  housed  in  a light  tight  case  which  also  contains  discriminator 
and  pulse  shaping  electronics.  The  output  from  the  detector  consists  of  a train 
of  identical  pulses  whose  mean  repetition  rate  is  proportional  to  the  intensity 
of  the  light  incident  on  the  photo-cathode. 

The  third  part  of  the  optical  system  is  the  monitor  arm;  this  consists  of 
a similar,  but  fixed,  arrangement  to  the  moving  arm.  It  is  positioned  at  a 
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scattering  angle  of  approximately  60  degrees,  but  on  the  opposite  side  of  the 
illuminating  laser  beam  to  the  moving  arm.  Together  with  readings  of  the  laser 
power  meter,  it  provides  the  means  of  correcting  for  the  fluctuations  in  the 
number  of  the  particles  passing  through  the  scattering  volume. 

For  convenience  the  output  pulses  from  the  photomultipliers  were 
counted  by  the  Malvern  digital  autocorrelator,  type  K7023.  This  had  the  advantage 
that  it  simplified  the  experimental  procedure.  For  each  angle,  the  total  number 
of  photon-detections  could  be  counted  over  a preset  period  and  automatically 
recorded  on  paper  tape  for  subsequent  analysis. 

The  particle  generation  and  filtering  arrangement  is  shown  in  Fig  13.  It 
consists  of  a commercial  oil-mist  generator  (Norgren  micro-fog),  intended  for 
use  in  lubricating  systems,  and  an  optional  5 litre  settling  chamber.  This 
could  be  used  to  remove  large  particles , formed  by  agglomeration,  and  to  smooth 
out  fluctuations  in  the  number  of  particles  passing  through  the  sample  volume. 

The  mean  rate  at  which  the  particles  were  produced  was  controlled  by  adjusting 
the  drip-rate  and  hence  the  volume  of  the  oil  supplied  to  the  atomising  venturi 
section. 

The  seeds  were  delivered  to  the  scattering  region  at  a pressure  slightly 
above  atmospheric  and  launched  into  the  atmosphere  through  a narrow  glass  tube, 
approximately  1 mm  in  diameter,  as  shown  in  Fig  II.  In  order  to  avoid  contamina- 
tion of  the  optics  and  to  prevent  the  accumulation  of  oil  droplets  in  the  labora- 
tory, the  particles  were  drawn  away  from  the  scattering  region  by  an  exhauster; 
this  also  provided  a moving  sheath  of  air  around  the  scattering  region  which 
limited  the  ingress  of  dust  particles  from  the  laboratory. 

The  experiment  was  carried  out  by  recording  the  total  number  of  photon- 
detections  (counts)  from  the  detectors,  in  some  convenient  time  (10  seconds)  at 
each  scattering  angle  0 . This  was  repeated  a number  of  times,  at  each  angle, 
to  enable  the  elimination  of  abnormally  high  counts  to  be  carried  out;  these 
anomalies  could  arise  from  example  when  a large  agglomerated  particle  passed 
through  the  scattering  volume. 

A complete  set  of  experiments  was  performed  over  the  entire  range  of 
scattering  angles  0 ( 1 0- 1 60  degrees)  in  5 degree  steps  and  a comparative  set 

of  experiments  was  carried  out  for  the  cases  where  the  particles  were  absent  and 
the  laser  beam  switched  off.  This  permitted  corrections  to  be  made  to  the 
scattering  results  for  the  effect  of  background  radiation. 
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3.3  Extinction  measurements 

This  experiment  was  carried  out  in  an  attempt  to  verify  the  results  for 
particle  size  obtained  from  the  angular  scattering  measurements  described  above. 
However,  the  interpretation  of  the  measurement  depends  on  the  assumption  that  the 
micro-fog  generator  produces  particles  with  a narrow  size  distribution  (ze  a near- 
monodispersion). In  addition  the  technique  requires  the  experimentally  derived 
points  (log^fg/f  plotted  against  v ) to  be  matched  to  the  curve  Qext-(a)  • 
However,  for  the  range  of  a over  which  the  experiment  can  easily  be  carried 
out  this  method  is  subject  to  large  errors.  In  addition  there  is  the  conflicting 
requirement  that  the  experimental  range  of  wavelengths  and  hence  the  range  of 
values  of  a should  be  limited,  or  the  dispersion  of  the  oil  (change  in  refract- 
ive index  with  wavelength)  will  cause  a significant  change  in  the  curve 

Qext(°0  5 see  Fig  1. 

The  extinction  experiment  is  designed  to  measure  the  extinction  coefficient 
t of  a suspension  of  oil  droplets  (produced  by  the  oil-mist  generator)  at  a number 
of  different  wavelengths  and  simultaneously  to  measure  the  relative  number  of 
particles  present  in  the  sample  volume.  This  was  carried  out  with  the  arrangement 
shown  in  Fig  ]4. 

Two  laser  beams,  at  different  wavelengths,  pass  coaxially  through  the  drop- 
let suspension;  the  attenuation  of  one  of  the  beams  is  used  to  monitor  the 
relative  number  of  droplets  per  unit  volume  N , whilst  the  second  beam  gave  the 
relative  attenuation  at  a different  wavelength  for  the  same  droplet  density.  The 
optical  system  thus  has  three  functions:  to  provide,  simultaneously,  a pair  of 
coaxial  collimated  beams  at  different  wavelengths;  to  provide  an  enclosure  for 
the  droplet  suspension;  and  to  provide  a receiving  system  which  accepts  only  the 
unscattered  light. 

The  first  function  is  fulfilled  by  combining  the  outputs  of  an  argon-ion 
laser  and  a helium-neon  laser  at  the  interface  of  a beamsplitter  B as  shown 
in  Fig  14.  The  resulting  coaxial  beam  is  expanded  and  collimated  by  the  arrange- 
ment Lj,  . The  diameter  of  the  collimated  beam  was  approximately  10  mn.  The 
power  of  each  laser  beam  was  monitored  by  detectors  Dj  and  via  the  dis- 

persing prism  Pj  . The  wavelengths  available  for  use  were  632.8  from  the 
He-Ne  laser  and  457,  488  and  514.5  nm  from  the  argon-ion  laser. 

The  second  function  of  the  optical  system,  the  enclosure,  was  provided  by 
a 1.5m  tube  made  from  resin-bonded  paper.  It  was  150  mm  in  diameter  and  fitted 
with  perspex  ends,  holding  removable  glass  windows.  Oil  mist  from  the  generator. 
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with  or  without  the  settling  chamber  in  circuit,  was  introduced  into  the  tube  at 

a tangent  to  the  wall,  in  order  to  provide  a homogeneous  distribution  of  droplets 

within  the  sample  volume.  The  mist  was  exhausted  from  the  opposite  end  of  the 

tube  and  expelled  from  the  laboratory  to  prevent  the  accumulation  of  droplets 

in  the  transmission  paths  external  to  the  test  cell.  To  provide  an  overall 
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attenuation  of  the  transmitted  beam  of  80%,  as  recommended  by  Hodkmson  , it 
was  necessary  to  reflect  the  beams  back  through  the  optical  cell,  providing  a 
total  path  length  of  z = 3.0  m. 

The  coaxial  beams  emerging  from  the  cell  were  split  into  two  parts  by  the 
Amici  prism  and  passed  directly  to  the  receiving  optics,  which  consisted  of 

a focusing  lens,  aperture  and  detector;  these  components  were  selected  to  fulfil 

the  third  function,  of  angular  separation  of  the  transmitted  and  scattered  light. 
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Hodkmson  claims  that  if  the  angle  subtended  by  an  aperture  at  its  associated 
lens  is  less  than  one-tenth  of  the  angle  of  the  first  minimum  of  the  Fraunhofer 
diffraction  pattern  resulting  from  the  largest  particle  in  the  scattering  volume, 
then  the  angular  resolution  of  the  receiving  system  is  adequate.  The  combinations 
of  lens  and  aperture  - A^  and  - A^  used  in  this  study  gave  sufficient 
resolution  with  this  criterion  for  particles  up  to  7 microns  in  diameter  (for  a 
wavelength  of  500  nm) . Note  that  providing  the  criterion  is  met  the  beams  do  not 
have  to  be  accurately  collimated. 

The  detectors  are  large  area  silicon  photodiodes,  type  BPY  13A; 

operated  in  the  reverse  bias  mode,  these  detectors  have  a linear  response  to  light 
intensity  over  the  range  of  interest. 

The  first  part  of  the  experiment  consists  of  selecting  the  appropriate 
wavelength  of  the  argon-ion  laser,  aligning  the  optics  and  measuring  the  currents 
through  the  photo-detectors;  this  corresponds,  with  the  cell  empty,  to  the 
measurement  of  the  incident  intensity  f^  for  both  wavelengths.  The  cell  is 
then  gradually  filled  with  oil  droplets  until  equilibrium  is  reached,  when  the 
number  of  particles  per  unit  volume  is  constant;  at  this  point  the  transmitted 
intensities  f are  measured  for  both  wavelengths,  simultaneously;  in  addition 
a note  is  made  of  the  monitor  diode  currents. 

The  experiment  is  repeated  for  different  combinations  of  wavelength  and 
filtering  (of  the  particles).  The  data  reduction  consists  of  subtracting  the 
dark  currents  from  the  detector  values,  standardising  the  value  of  f against 
fluctuations  in  laser  power  and  eliminating  the  particle  number  density  by 
deriving  a relative  extinction  coefficient  for  each  wavelength,  compared  to  the 
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He-Ne  extinction.  The  points  corresponding  to  the  relative  value  of  log^fg/f 
plotted  against  wave  number  (v)  are  then  fitted  to  the  plot  of  Q (a)  obtained 

from  Penndorf's  data  for  monodispersions  with  a refractive  index  of  1.486. 

4 RESULTS  AND  DISCUSSION 

4. 1 Refractive  index  measurements 

The  absorption  of  the  laser  beam  passing  through  78  mm  of  the  OM-11  light 
mineral  oil,  at  approximately  22°C,  was  found  to  be  approximately  2.5%  at  a wave- 
length of  488  nm.  The  value  of  the  absorption  coefficient,  defined  in 

__  J 

equation  (19),  was  therefore  0.32  m , which  leads  to  a value  of  1.26  x 10  for 
the  imaginary  part  of  the  refractive  index,  n^  . 

The  value  of  the  real  part  of  the  refractive  index,  n^  , was  measured  at 
several  wavelengths  and  calculated  for  the  remaining  laser  wavelengths.  The 
results  are  summarised  in  Table  1. 

Table  1 


Real  part  of  refractive  index  n^ 

1.475* 

1.477 

1.479 

1.481* 

1.483* 

1.486* 

1.488 

Wavelength  nm 

632.8 

577 

546 

514.5 

488.0 

457 

436 

* Indicates  estimated  value  (using  Cauchy  relationship) 

It  can  be  seen  from  the  table  that  the  value  of  the  refractive  index  changes  by 
less  than  1%  over  the  wavelength  range  436-632.8  nm,  so  that  the  calculated 
curve  for  (^^(ci)  for  t^e  value  1.486  is  adequate  for  the  interpretation  of 
the  results  from  the  extinction  experiment,  since  other  errors,  for  example  from 
curve  fitting,  are  more  significant. 

4.2  Angular  scattering  intensity  measurements 

Figs  15  to  16  show  the  curves  I j (6 ) , 12(6)  obtained  from  the  experiment 
described  above.  As  can  be  seen  the  curves  show  none  of  the  detailed  structure 
predicted  by  the  Mie  theory.  This  suggests  that  the  particle  population  is  a 
polydispersion,  although  the  degree  of  dispersity  required  to  produce  such  smooth 
curves  is  not  known.  Some  insight  may  be  gained  by  inspection  of  the  curves 
drawn  in  Figs  3 to  6,  which  show  that  small  variations  in  the  particle  size  para- 
meter a produce  large  changes  in  the  structure  of  the  curves. 

It  can  be  seen  from  the  measured  distribution  that  a laser  anemometer 
collecting  light  in  the  near  back-scatter  direction  6 — 160  degrees  receives  far 
less  signal  than  in  the  corresponding  forward  scatter  direction.  The  minimum 
intensity  of  the  scattered  light  from  the  particles  used  at  RAE  appears  to  lie 
at  around  120  degrees.  055 
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The  angular  dependence  of  the  scattered  light  intensity  for  unfiltered 
particles  (Fig  16)  shows  relatively  greater  forward  scattering  than  in  the  case 
of  the  filtered  particles  (Fig  15).  This  is  consistent  with  the  presence  of  larger 
particles  in  the  scattering  volume. 

The  derived  Sloan-type  plots  for  the  results  obtained  with  the  settling 
chamber  are  drawn  in  Fig  17  for  the  two  polarisation  states.  The  curves  have  clear 
first  maxima  at  angles  of  38  degrees  and  44.6  degrees  which  yield  from  equation  (II) 
values  of  0.48  and  0.41  microns  for  the  mean  diameter  estimates.  On  the  assump- 
tion that  the  particles  are  a true  monodispersion,  the  Mie  theory  based  correction 
for  the  refractive  index  and  particle  size  (Fig  7)  would  give  corresponding 
diameters  of  0.51  and  0.44  microns.  These  values  are  consistent  with  experience 
gained  in  the  laser  anemometry  experiments  in  rapidly  accelerating  flows  referred 

to  in  the  introduction,  and  also  agree  with  an  independent  investigation  of 
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particle  size  carried  out  by  Gregory 
4.3  Extinction  measurements 

Fig  18  shows  the  limited  number  of  points,  derived  from  measurements  of 
the  relative  extinction,  fitted  to  the  curve  of  Q t(a)  • It  clearly  shows  the 
difficulty  encountered  in  interpreting  the  measurements;  however  the  correspondence 
of  the  abscissae  indicates  that  the  particle  diameter  is  approximately  0.7!  microns 
for  the  filtered  particles.  Measurements  made  on  the  unfiltered  particles  indi- 
cated that  they  had  a rather  larger  diameter  (0.74  microns). 

5 CONCLUSIONS 

The  results  of  the  experiments  confirm  that  the  seeding  particles  used  at 
present  for  wind-tunnel  laser  anemometry  at  RAE  are  suitable  for  most  applications 
of  interest.  The  results  obtained  from  the  experiments  show  that  the  particles 
are  well  within  the  diameter  range  0.1  to  1.0  microns,  which  accords  with 
experience  gained  in  wind-tunnel  experiments  on  strongly  accelerating  flows.  The 
values  around  0.5  micron  for  the  mean  diameters  obtained  from  the  Sloan  plots  are 
more  reliable  in  view  of  the  difficulty  of  interpreting  the  other  forms  of  data. 

The  experiments  have  also  shown  that  large  particles  can  agglomerate  in  the 
seeding  system  and  that  these  must  be  filtered  out  if  misleading  estimates  of 
velocity  in  highly  accelerating  flows  are  to  be  avoided.  It  has  also  been  shown 
from  the  shape  of  the  angular  intensity  distribution  of  the  scattered  light  that, 
for  anemometry  experiments  using  this  type  of  seeding,  scattering  angles  around 
120  degrees  should  be  avoided. 
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The  overall  conclusion  reached  in  this  study  is  that  the  accurate  determina- 
tion of  the  size  of  the  particles  whose  diameter  is  similar  to  the  incident  wave- 
length is  very  difficult  using  light  scattering  techniques.  Although  more 
detailed  size  experiments  on  the  oil  droplets  are  probably  not  necessary  for 
existing  anemometry  applications  at  RAE,  an  experiment  of  the  Millikan  type 
where  measurements  can  be  carried  out  on  a single  particle  would  be  very 
instructive.  Methods  based  on  Fraunhofer  diffraction  at  ultra-violet  wavelengths 
could  also  possibly  be  developed  to  yield  estimates  of  the  size  distribution. 
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Fig  8 Measurement  of  the  imaginary  part  of  the 
refractive  index  (m) 


Reflected  ray 


Hollow  prism 


Plan  view 
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Fig  10  Experimental  arrangement  for  measurement  of  real  part 
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Fig  1 1 Schematic  diagram  of  angular  scattering  experiment 
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Fig  14  Optical  arrangement  of  extinction  experiment 


